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Abstract 


The  iocoMplete  Dlrlchlet’a  Multiple  Iptegral  Is  defined  as 

a 

► 

2 


•••  j f(ti  + t2  + •••  + t^)  t^  t^^  •••  t^  dti  dt2 


dt. 


where  R - {(t^,  t2*  •••.  t^) : i 0,  | < 1,  t^^j^  <_  a^  ^ t^;  1 - 1,  2,  •••,  n - 1} 


It  Is  shown  that  this  Integral  Is  reducible  to  a single  Integral  Involving  products 
of  inconplete  beta  fuctlons  as  a multiplicative  constant.  The  results  are 
generalised  and  extended  to  other  aniltlple  Integrals.  Finally,  as  an  application 
a class  of  multivariate  type  distribution  defined  by  a functional  form  Is  par- 
tially Investigated. 


Introduction 


If  £(•)  is  a continuous  function  of  some  real  non-negative  arguawnt,  and 
*t  n,  then  the  Dlrlchlet's  multiple  Integral 

I I - I -1  • -2 


> 0,  !■  2| 


t,+t-+***+t  <1 
1 £ n— 


a^-1  <»2“1  o “1 

f(ti  + t,  + •••  + tjj)  tj^  t2  •••  t^"  dt  dt2  •••  dt^ 

(1) 


may  be  expressed  as  the  single  Integral  [5],  [18] 


r(o 

r(^ 


- + a,  + •••  + o ) i 


(a,+02+‘**+aj^)-l 


dr 


*1  "2  "n'  0 

Alternatively,  we  may  express  the  multiplicative  quotient  Involving  gan 
tlons  in  terms  of  products  of  beta  functions  and  write 


func- 


Iq  ■ 8(02*  Oj^)B(aj,  O2  + oij^)  •••  B(a^,  + o^_2  + •••  + O2  + Oj^) 

i: 


(o,+a-+**»+a_)-l 
f(T)T  "dr 


(2) 


%rhere  In  general 


An  extension  of  (1)  was  proposed  by  Slvasllan  [16]  where  in  (1)  the  argument 

of  the  function  f(<)  was  replaced  by  various  partial  sums  of  the  variables  tj^, 

^2*  ***’  ^n’  extensions  In  the  Integrand  form  of  (1)  were  proposed  by 

Slvazllan  [17]  and  Klamkln  [13).  Klamkln  also  suggested  an  Interesting 

extension  of  (I)  in  yet  another  direction,  namely  by  changing  the  region  of 

r n 

integration  to  {(t^,  tj,  ....  t^^):  1 0,  J t^  < t,  J t^  < t'},  t,  t'  > 0. 

1 r+1 

The  multiple  integral  le  then  reducible  to  a double  integral  again  Involving 
quotients  of  gamma  functions  as  Its  multiplicative  constant. 

In  the  present  work,  we  propose  still  another  region  R defined  for 
0 < a^  £ •*,  1*1,  2,  . . . , n as  the  simplex 


1 


¥■  - 


'.t 


«■  A 

R » t2 t^):  > 0,  I < 1,  1 I 

1-1,  2 n - 1}  (3) 


and  define  the  lnco«a>lete  Dlrlchlet’s  multiple  Integral  as 


Hf  ®‘l~^ 

•••  j f(tj  + t2  + •*•  + t_)  t,  t. 


a -1 
^ n 


* dt,  dt,  ...  dt  (4) 

1.  £.  n 

This  integral  is  not  reducible  to  any  of  the  previous  forms  investigated. 


However,  in  the  very  special  case  when  a^^  - ®,  the  Liouville's  extension  of 
Dirichlet's  integral,  that  is  form  (1)  is  recovered.  In  general,  we  shall  show 
that  the  proposed  multiple  integral  (4)  has  the  same  format  as  (2)  except 
that  the  multiplicative  constants  outside  the  single  integral  are  products  of 
incomplete  beta  functions.  We  shall  show  that  - kl^,  where  k is  the  product 
of  (n  - 1)  incomplete  beta  function  ratio. 

2.  The  Incomplete  Dirichlet's  Multiple  Integral 

2 

To  reduce  (4),  we  first  make  the  change  in  variable  1 “ 1»  2,  ...i 

n.  The  Jacobian  of  the  transformation  is  clearly  z”xjX2 

A 

region  of  integration  R is  defined  by  the  set  of  points 

n 


The  new 


R = {(x^,  X2 Xjj):  Xj  > 0,  I *1  1 *1+1  i + X2  + . . . + xj 

for  1 - 1,  2 - 1} 


Hence 


»» • n I '<’!  * 


2 2 ^Oj^-l  2o2-1 

X2  ^ . . . ^ *n^  ^1  ^2  • • • 

dx,  dx-  . . . dx 
12  n 


2o  -1 
n 


He  now  utilise  a second  change  in  variable  which  is  the  elegant  transformation 
into  generalised  n-diawnaional  spherical  coordinates  proposed  by  Schlafll  in 
1855  [4],  also  by  Clare  in  1881  [51,  and  used  often  in  statistics  (e.g.  see 
[12]).  Let  then 
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X , ■ r cos  6 , sin  0 - 

n-1  n-1  n-2 

X ■ r sin  0 , 

n n-1 

The  Jacobian  of  the  transformation  Is 

r"  ^(cos  0 ^(cos  0 ^•••(cos  0.) 

n— i n— z i. 

To  determine  the  new  region  of  Integration  R,  ve  note  that  the  ball 

n 2 

0 £ X £ 1 Is  mapped  Into  the  set  of  points  0 _<  r £ 1.  Also,  for  1*1,  2, 
1 ^ 

n - 1,  the  set  of  points  defined  by 

0 £ *^1  *2  ^ *’*  *1 
Is  mapped  Into  the  set  of  points 

0 < r cos  0 , cos  0 - • • • cos  0 , . _ cos  9 . . , sin  0 . 

— n-1  n-Z  1+2  1+1  1 

< /T"r  cos  0 , cos  0 o CCS  0..,  cos  0. 

— 1 n-1  n-z  1+1  1 


or 

or 


0 £ tan  0j^  £ 


0 £ 0j^  £ tan 


It  thus  follows  that  R Is  mapped  Into  R where 

R - {(r,  0,,  0_,  ....  0 0 < r < 1,  0 < 0.  < tan 

i i n— 1 — ■“  — X — X 

for  1 ■ 1,  2,  ....  n - 1} 
The  expression  for  i.  can  then  be  written  as 


t - 2 


2 20j^-l 

f(r  ) (r  cos  0 , cos  0 -•••  cos  0,  cos  0.) 

n-1  n-2  2 1 


(r  cos  6 - cos  0 -•••  cos  0_  sin  0.) 

n— 1 n-1  c X 

-1 

•••  (r  cos  0 - sin  0 ,)  (r  sin  0 ,) 

n“i  n*i 

Ir"”^(cos  ejj_j^)“”^(cos  ^***  **®1  **®2' 


I \ 
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r 

- 2 I 


fa 


• *d0  I 

"I 

( 

I 

« 


tan 


7 “1  2a  -1  2a- -1 

t j (sin  0j^)  (cos  0j^)  ^ d0j^] 

tan 

f 2 203-1  2(o,+a  )-l 

Jo  ®2^  **®2^ 


tan 


tan 


■ ““n-l”^  2(a,+o-+. . ,+a  ,)-l 

(sin  0^,)  " ^ (cos  0 ^ ^ "-2 


‘*V2l 


f ^ ^ ^®n~^  2(a^+a-+. . .+0  )-l 

' Jo  Vl>  ” 


Making  the  change  in  variable  x » and  ■ sin^  0^,  1 - 1,  2,  ..., 


n - 1,  we  obtain 


fl+ai  a^-l  ,1+a  o.-l  a,+a,-l 

n “ f 1^  “1  - “l>  «*“il  [ ■ “2^  ^ ^ **“2^ 


n-1 

rl+a_  ,0-1 


- ' Jo  - VP  "■*  ° 

if  t(T)T°^^^* 


Vl'^n-2'^**-'*^2'^r^ 


‘‘ViJ 


Using  the  definition  of  the  Incomplete  beta  function 


B,(a,  8)  - y®"^  (1  - y)®-l  dy 

^0 


we  obtain  the  following  expression  for  I 


“ «j  <“2-  V”*  a,  <°3’  “2  * “P> 


“Vl  Vl  * 

^+*n-l 

V2-"  ••*  +“l>J 


f(T)T  ^ ^ dr 

^0 


4' 


Finally,  If  b)  “ b)/B(a,  b)  is  the  Incomplete  beta  function  ratio, 


then  clearly: 

n-1 

I /I  - n I 
" " 1-1 

1+a. 


(«!,  + o^_2  + 


which  is  independent  of  the  function  f(*)*  In  the  special  case  when  ■ <*> 
for  all  i-  1,  2,  ...,  n-  1,  the  value  of  the  previous  expression  equals  1. 


itim 


where  T is  the  region  of  integration  defined  by  the  set  of  points 


V ® 

T - {(v.,  V , ....  V ) : V > 0,  ^ (-q^)  £ t, 

11  ^1+1 


r ''j 

— ®i  L ^0  ^ ;i»l,  2,  ...,n} 

1 j 


under  the  assumptions  that  > 0,  > 0,  > 0 for  1 ■ 1,  2,  ....  n. 

Vf  ®1 

making  the  change  In  variable  (— ) - u. , 1 - 1,  2,  ....  n,  we  obtain 

1 ^ 

“l  “2 

q,  q 


K„(t) 


■‘b  6 " f f I “"l  * "2  •••  * 

_1_1  -1-1  _n 

6,  B,  B 

12  n , , , 

u,  u„  •••  u du,du„  •••  du 
1 2 n 1 2 n 


which  is  of  the  same  form  as  (6).  Hence,  using  (7)  we  obtain 

a 


“l  “2 

q q * 

1 2 


B. 


[B 


a a 

<b!- b:>"»  .2 


l+ai 


'2  ^1 


1+a- 


(!l  ^2 

S*  ®2 


••IB  an_i 


“ ® -1  o o a, 

/ o n*x  1 , X\ 1 

,t  6j  6j  B„ 


Jq 


f(T)T 


Here, 


dr 


(9) 


4. 


Extensions 


0,-1  02*1 


o -1 

t " dt,  dt,  •••  dt 
n 12  n 


We  now  consider  for  r - 1,  2,  n and  0 s £ n - r,  the  Integral 

' I I ’1‘  I +1  + *•*  + Vs>'^/  *^2 

K 

Where  R is  the  sane  region  as  described  in  (3).  The  function  f(*)  and  the 
parameters  o^,  02»  ..•>  o^  are  as  defined  before.  It  Is  possible  to  reduce 
- and  J - to  single  Integrals  although  a simple  reduction  does  not 
appear  to  exist  for  other  values  of  r and  s 
a.  The  Integral  ^_2 

We  first  show  that  a double  integral  of  the  form 

H(a,  b)  » (q  (q  f[T(l  - v)AHt(1  - v)]“~^  (tv)®  ^ t dx  dv 

where  0 _<  a,b  £ 1;  A > 0;  o,E  > 0 Is  reducible  to  single  Integrals.  Making 
the  change  in  variable  t(1  - v)  “6  and  (1  - v)  ■ u,  we  get 

rau 


b)  - /i_b  /o 
Interchanging  the  order  of  Integration  yields 


f(A0)e“‘*’®‘^  (1  - u)®~^  u®*^  de  du 


H(a,  b)  - jQ  f(A6)e“^®’^  (1  - u)®"^  u’®“^  du  d6 


- ’’(1  - u)®"V®"^ 


du  d6 


Making  the  change  in  variable  t ■ (1  - u)/u  we  obtain  finally 
(•a 


Ka.  b)  - f /q  f(A6)0“^®“^(|  - D®  d0 


■ B *’^(A0)0®‘^®"^I(f  - D®  - (j-rt)**]dO 


Now 


-'l.n-2  “ ^ ^ ^^2^  •••  •••  C 


a -1 


dt,  dt,  •••  dt 
12  n 


To  evaluate  5 . we  perform  the  same  sequence  of  variable  transformation 

X f 


as  In  the  evaluation  of  '1^,  and  obtain 


“^l.n-2  ‘ t®  ai  <“2*  “l>JI®  *2 


l4a. 


l+a* 


(“3*  “2  “l>^  •••  a„-2 

K-V  V2  + “n-3  + •••  + “i)l 


n-1 


/1+a  , [1 

I I 

h h 


f[T(l  - u^_i)][t(1  - u^_p] 


a,+o„+***+a  --1 

1 £ n~i 

a -1 

(t  u 1 ) " T dx  du 


n-1'  ‘ “n-1 

This  double  integral  Is  of  the  same  form  as  H(a,  b)  with  A - 1,  a > 1, 

b * a ,/(l  + a 0-0,  + a,  + •••  + o ,,  0 - o^.  Hence 

n-i  n— 1 1 z n-1  n 


«n-2 


"^l.n-^  “ f®ai  <®2’  a2  ^®3»  <*2  ■*■  “l^^  *”  ^® 

TFi^  l+«2  l+*n-2 

<Vl*  “n-2  “n-3  + -•  + “i)l 


1 

l^lo  f(e)e 

n 

1 

fl+a 


*^T —X  « V* 

^ “ (i  - 1)  “ de 


, /A‘*‘a_  , o,+o,+»**+o  -1  , » 

- ^ Jo  f(e)e  ^ " [(j  - 1)  “ - (a^_i)  "ide 


and  In  final  form 


’l,n-2 


I*  «i  ^*2*  “1^11®  a2  ^“3’  ®2  “1^^  ***  ^® 


l-fa2 


•n-2 


TF* 


n-2 


‘Vr  “d-2  * V3  + V3  * •••  * “l” 


9 • 

^wsew'^rTaiTigs''— B 


u/;  - 

a 1 

(a  " rl+a  , 

+ "-^£<9) 


a,+o-+***+a  ,-l  a 

I ^ ^ (1  - e)  " de 


rl+a  , o-+a_-«-**+a  -1 

j f(e)0  ^ ^ “ de} 


b.  The  Integral  q 

We  consider  here  a double  integral  of  the  forie 

K(a,  b)  - /q  /q  f(TvA)[T(l  - v)]“"^(tv)®"^  t dr  dv 
where  0 ± a^^b  <.1;  A > 0;  a,  0 > 0.  Making  the  change  In  variable  tv  ■ 0 


and  V " u,  we  get 


K(a,  b)  - /q  /q  f (A0)e“'*‘®"^(l  - u)“”^  u ® ^ d0  du 

- /q  f(A0)0®'^®‘^  /g  (1  - u)“"^u"“"^du  de 

a 

- i f(A0)e“'*'®"^[-(^  - i)°  + (f  - 1)  ]de 

“ " o^i  ■ Jo  ‘*®  ■*■  o ^0  f ■ ®>“  **® 

Now 

nf  0,-1  0,-1  0-1 

•••  I £„  -“i  "'2  ••• 

R 

Again,  we  perform  the  same  sequence  of  variable  transformation  as  in  the 


evaluation  of  2.  and  we  obtain 
n 


*1 

TOT 


(02.  o^)][B 

TRKT" 


<V  «2  + V>  ••• 


*n-l 


f A+*n-l 

<®n-l*  ®n-2  ®n-3  **'  * ®1>^  h >0  ^<"Vl^ 

o,-fo,+«*«+o  ,-l  0-1 

[t(1  - u _,)1  ^ ^ ^ (tu_,)  t dt  du^, 


This  double  integral  Is  of  the  sane  form  as  K(a,  b)  with  A - 1,  a • 1 


-a  ,/(l+a  +a,  +***+a^.  and  8 - o . Hence 

n— 1 n~i.  X * n~x  n 


'•n.O  • .1  <“2-  «2  '“3>  “2  * “P'  •••  '*  ,^,.2 


1+ai 


TTaJ 


1+a, 


n-2 


^Vr  “n-2  “n-3  "* 


I [ f 

a,  + o_  + •••  + a , ‘ '0 
1 Z n— 1 


n-1 

1+a  , a -1  o +a-+***+o  . 

f(0)e  " (1-0)  d0 


o  +o»+* • *+a  , n-1 

1 1 X n.i  , a +a,+*“+a  -I 

f(0)0  ^ 2 


d0] 


5.  A Still  Further  Extension. 

Another  multiple  integral  tihich  is  reducible  to  a single  Integral  is  the 
following 


II  •••  I + ^2  "^  ***  "^  ^n^  ‘ 

R ^ 


a,-l  a--l  a -1 

^ ...  t ** 

2 'n  . 


0 0 0 
(ti  + t2>  ^(tj^  + tj  + tj)  ^ (tj^  + t2  + •••  ^ 

dti  dt,  dt 
■*■2  n 

where  8^  1 0,  (1*2,  3,  •••,  n - 1).  The  reduction  process  Is  slsillar  to  that 
used  for  (4)  end  will  not  be  repeated  here. 


11 
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6.  Applications  to  Statistics 


Statistical  distributions  of  the  multivariate  type  classified  by  functional 


form  have  received  little  attention.  For  univariate  classes,  the  Pearson's 


type  (e.g.  see  [6])  are  well  knom.  For  the  multivariate  case.  Lord  [14]  and 


Box  and  Hunter  [3]  consider  the  radical  or  spherical  distributions  where  the 


joint  probability  density  functions  of  the  sequence  of  random  variables 


where  A is  a constant.  Kelker  [11]  studies  linear  transformations  of  spherically 


distributed  variables,  while  Baldessarl  [2]  considers  density  functions  of  the 


- *n>  J 

n n 

2 2 

where  I (x.  ~ <.  r • Finally,  the  so-called  elllptlcally  symsmtric 

J-1  J J 

distributions  are  studied  by  He Grew  and  Wagner  [IS]. 

Power  transformations  of  the  type  X ■ have  been  used  In  Dlrlchlet's 


distributions  (see  e.g.  Johnston  and  Kots  [9])  and  In  the  form  X ■ aY  - (1  - Y) 


are  known  as  the  TMcey  lambda  transformation  [8] [10].  Generalisatlona  to 


the  standard  Dlrlchlet  distribution  were  considered  by  Johnston  and  Kots  [9] 


We  consider  presently  some  generalisations  of  the  previously  investigated 


classes  of  statistical  distributions.  In  particular,  we  study  soma  properties 


of  a class  of  multivariate  distributions  defined  by  a functional  form  in 


which  Che  variables  aay  be  subject  to  a power  transforaatlon. 

Let  then  A be  a positive  constant  and  f(x)  be  a continuous  nonnegacive 
function  defined  everywhere  in  the  interval  0 £ x < • with  the  property 
that  for  Y-  >0,  i - 1,  2,  ....  n,  the  ioproper  integral 


exists  and  converges  absolutely 


X } be  an  n-dinensional  non-negative  random  variable 


whose  P.D.F.  ^ 


over  Che  region  { (x. 


n the  quantities  a.  > 0,  g.  > 0 


q.  > 0 be  uniquely  defined 


a.  The  Constant  A 


The  quantity  A is  given  by  Che  normalising  conditions 


The  multiple  integral  is  of  Che 


b.  Mixed  Moaenta 


The  alxed  aoBents  about  Che  origin  if  they  exist,  are  defined  for 
1,  2,  (1  ■ 1,  2,  n)  by  the  expression 


c.  Marginal  Distributions 


The  Marginal  PDF  4 


By  a suitable  change  in  the  order  of  integration,  setting  y ■ t (x^/q  ) 


and  using  the  convolution  theorea,  one  can  verify  that 


siallar  fashion,  it  is  possibla  to  shov  that  cha  joint  Marginal  P.D.P.  of 
(Xj^,  X2*  ••••  X|^}  (k  £ n)  is  given  by 


Cl  a 

“fc+l  „“k+2 


♦x  X X (*i*  *2’  •*•» 
1*  2“--’  k ^ ^ 


\:r  vr  v “2-^ 


®fcfl®k+2 


1— *1  *2  •••  *k 
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fi  A ^JStk  X “fc+2  . . . “n  , 

r»  V "o  » “«,  fl *S ••  + 1 

fl(^)  + (^)  + “•  (^)  +r]  d 

o \ ^2  \ 

The  conditional  distributions  can  than  readily  be  obtained.  We  presently 
extablish  the  following  two  theoreaw 
d.  Theorem  1 

The  (n  - 1)  random  variables  {Y, , Y-,  ....  Y ,}  where 

i.  £ n~i 


’*  " *1  ®1  *,  ®2  » 'n 

(^)  + (^)  + — + (^) 

12  n 


• l*lt  2f  ...f  n*l 


have  a Joint  probability  density  function  which  is  of  the  Dirichlet  type 
inside  the  simplex 

I ■ Htj.  Y2.  •••t  y„_i)  : yj  + y2  + •••  + < i.  y^  >.  o. 


X2  1 0.  •••.  y„_i  t 


and  which  is  xero  elsewhere. 


Proof; 

*1  “*  *. 

’«  ■ <5T>  ♦ •••  ♦ <q^> 

1 ^n 

Then  (-jp)  - Y^Y^^,  i - 1,  2 - I 

- I'l  - »2 Vl’*. 


This  particular  tranafonoation  aaps  the  region  {(x^,  Xj,  ....  x^^)  : 0 < x^  < •, 

i “ 1.  2 n}  onto  the  region  C - y^,  ....  y^)  : Xj  + V2  + •**  + y^_j  1 1, 

i 0 for  1*  1,  2,  n-1,  O^y^  <•}.  The  corresponding  Inverse 

functions  for  this  one  to  one  transforaatlon  Is 

1/B. 


2 f eeag  H“X 
l/6_ 


*n  ■ ’\.«1  - *1  - »2  - — Vl”.’  ° 

The  Jacobian  of  the  transformation  Is 


J - 


3(Xj^,  ^2 • •••»  Xjj) 


®(yj^»  X2*  •••»  y^) 


3x, 


3x, 


3x, 


3yi 


3x  3x 
n n 


3x 


3yi  3y<i 


3y. 


Now  for  1 ” 1,  2,  n,  we  have 

1 


3. 


- 1 


’*1 . , r‘ 

8^  <Vn>  ^ 


3x 


■^-1 


— 1 1 1 

■ IT  • #7  <V«>  »i 

n 1 


3x^ 


8y 


J 


- 0 for  J ^ 1,  n 


Also  for  1 ■>  1,  2,  ....n-1 


**«  ‘•n 

J®  ^ o 


- Vl>'n) 
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■ - . ■■ ■ ■ 


»s  q 1 

^ (1  - ri  - yj yn-l>f<l  “ - ^2 “ 

n Q 


The  value  of  tha  Jacobian  la  thus 

1 , 


\ S ®2  ‘‘n  1 ^ 

^ <V.>  • r <»2»n>  — FT  'VlJ-n* 

1 & n—1 


f-' 


r * ^2 vi>ynJ  ” 


0 y. 


-^n  -^n  * • • Vl> 


Noting  that  the  value  of  the  deternlnant  la  y^  , we  have 


— - 1 -i-_  1 

V2  •••  Jl  h 


Th  ^2  — Vi  y« 

n 

i-i 


1 , _L  1 

\-l  Jl  ®2  *'*®n 


- ^2 Vl>  " 

Hence,  the  joint  P.D.P.  of  (Y^,  ....  Y^}  la  for  (y^,  y^,  ....  y^^)  e 

-i-l  «2-l  a^^-1 


B.  «,-l 


a ,-l 


A (yiy„)  ' • <12'  (y2y„)  ' - Ci'  (Vi^n^ 


V ‘ yi  * T2 Vl^^nJ 


a,  a 
1 „ 2 


- A 


V •••  ^ 

»1®2  - ®n 


- 1 


■-  1 


n-1 

Vi 


...  . (1  . . 


“l  . “2  4.  “n 

sT^r*— 

'o-.'V 


--  1 


- Vi> 


and  is  zero  elsewhere.  Using  the  computed  value  of  the  constant  A In  (11)  we 
obtain  for  the  Joint  P.D.F.  of  {Y^,  ¥2.  Inside  the  simplex  ^ to  be 


“i  “9 

®1  h ®n 


!i.x 

®1  ®2 


--  1 


“1  “2 

h h 


a 

r(/) 

n 


(1  - Yi  - y,  - 


n-1 


•-  1 


- yn-l> 


and  to  be  zero  elsewhere. 


I 


e.  Theorem  2 


The  n random  variables  {Y. , Y, Y } where 

X 


(^> 
y . 1-H 

^ X.  h X,  ®2 


i+1 


1 ■ 1,  2 n - 1 


+ (-5^)  + •••  + (tt) 


X X ®2 

% ft  ^ 


X 6 
X n 


’n  ■ * ‘i;)  * - + 

are  Independently  distributed.  The  marginal  P.D.F.  of  Y^'s  are 


H-1 

y^^i+i 


'1 


♦y  ' 

*1  ^ 


^ ^ . . “n-l 

h h ***  ®i+i 

(1  + Yi)  ^ 

0 

“l  . ®2  . 4.  “n  . 

*n  ‘O-n)'. 


♦y  <Xn> 

n 


0 < < “» 

(1  - 1,  2. 

otherwise 


0 < y < c* 

'n 

otherwise 


where  the  A^'a  are  normalizing  constants. 
Proof: 


The  Joint  distribution  function  of  {Yj,  Yj,  ....  Y^^}  Is  given  by 

x<Xi  1 y^.  ^2  < Vj.  • • • t 1 y„>  - 

X ®2  X X »3  X h \ ®2 

P<(Tg)  . (-i^)  lyjK-^)  + 1 


X 


l"n-l  X i«n  X ®1 


(q->  lXn-1  ^ 

n " ^ 1-1  ^1  1-1  \ “ 
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Hence  the  P.D.F.  of  (1  ■ 1,  2,  n - 1)  is  for  0 < < • 
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